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The classical Lorenz lowest order system of three nonlinear or- 
dinary differential equations, capable of producing chaotic solutions, 
has been generalized by various authors in two main directions: (i) for 
number of equations larger than three (Curryl978) and (ii) for the case 
of complex variables and parameters. Problems of laser physics and 
geophysical fluid dynamics (baroclinic instability, geodynamic theory, 
etc. - see the references) can be related to this second aspect of gen- 
eralization. In this paper we study the asymptotic properties of some 
complex Lorenz systems, keeping in the mind the physical basis of the 
model mathematical equations. 



1 Introduction 

The behaviour of the classical Lorenz system of equations (Lorenz 1963) 



where a > 0, r > 0, b > 0, is well known (Shimizu and Morioka 1978, Yorke 
and Yorke 1979, Franceschini 1980, Sparrow 1982, McGuinnes 1983, Schmutz 
and Rueff 1984). Here X(t),Y(t), Z(t) are real functions of the time t and 
(a, r, b) are real parameters. The goal of this study is to ascertain to what 
extent some fundamental properties of the system are presented in its 
complex analogs. For comparison we choose the following properties of the 
system of equations ((H): 
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1. The fixed points of ((TJ 



(X,Y,Z) = (0,0,0), r>0, (2) 

(X = Y= ±\fti2,Z = r-l), r>l. (3) 

2. The conditions for stability of the second fixed fixed point in Q are 

a(a + b + 3) , i 

1 < r < r H = — -, b + 1 < a < oo. (4) 

o — b — 1 

At r > a chaos occurs. 

3. If cr — > oo then r# — > oo and the above fixed point remains stable for 
arbitrary r > 1. In this case 

Y(t)=X(t), 



d 2 Z „ „ 



(6-2/ + 2Z)— -26Z(/-Z)=0, (5) 



dt 2 y ' dt 

where / = r — 1, so that chaos is impossible. 

4. If 6 = 2cr and t — ■> oo, then 

d 2 X . . dX . . 1 q . . 

so that again chaos is impossible. The single second order autonomous 
equations in (JSJ and © are of known types (Lienard and Duffing re- 
spectively - Perko 1996). Their derivation can be seen below. 

The Lorenz system ([T} is a very simplified model of thermal convection in 
fluids (the Rayleigh-Benard problem, see Lorenz 1963 or Tritton 1988). Later 
on, problems of laser physics (see Ning and Haken, 1990 and the references 
therein) and of so called baroclinic instability of the geophysical flows (in the 
atmosphere or in the ocean) unexpectedly resulted in a system of the same 
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structure like (JIJ but in the complex domain (Fowler et al. 1983, Rauh et al. 
1996) 



dY 

= rX - aY - XZ, 

dt 

°— = -bZ + -(X*Y + XY*). (7) 
dt 2 

where i = r = r 1 + ir 2 ,a = 1 - ie,{X,Y) = (Xi,Yi) + i(X 2 ,Y 2 ), 

(a, b, e, Z) are real quantities and * denotes complex conjugate. Another 
similar system originated from geophysics will be introduced in section 4. 

In real variables the system (|7J) is equivalent to the following fifth order 
system of ordinary differential equations 



dX x 

~df 

dX 2 

dt 

dY\ 

~dt 
dY 2 



-aX x + aY u 
-aX 2 + aY 2 , 



r x X x -Y x - r 2 X 2 - eY 2 - X X Z, 
n X 2 -Y 2 - r 2 X x - eY x - X 2 Z, 



dt 

dZ 

_ = _Z + X 1 Y 1 +X 2 Y 2 . (8) 

In what follows we show that in some particular cases the above system 
simplifies and can be integrated even analytically. 



2 Particular cases of the system (Ej) 

2.1 Case a —> oo 

Let a — > oo, i.e., e = l/cr — > 0. Then the first equation of (JTj) implies 

Y(t)=X(t), (9) 



so that (0) reduces to 



dX 

— = (f-a-Z)X, (10) 
f = ~»Z + | A' |», (11) 



where | X \ 2 = XX* = X\ + X\ = A 2 > 0. In real variables the above equa- 
tions are equivalent to a three-dimensional system of nonlinear differential 
equations 

^ = (A - z)x, - ux 2 , 

at 
dX 

-J. = (f l -Z)X 3 + uX 1 , 

d ^ = -bZ + X 2 + X 2 , (12) 

where fi — r x — 1, v = e + r 2 . From the first two equations of ()12J) we obtain 
easily 

-T.A 2 = 2{h - Z)A 2 , (13) 
at 

Eliminating A 2 from (|13|) and (jllj) we obtain a single equation of known 
(Lienard) type (Perko 1996) 

d 2 7 d7 

— + (b - 2f x + 2Z) — - 2bZ{h -Z) = 0, (14) 

In addition by eliminating either X\ or X 2 from the first two equations of 
()12|) we obtain the second order equation of kind forced harmonic oscillator 



d 2 F dF 



0, (15) 



where F is equal to X\ or to X 2 . The equations (fT4*j) and ()15|) form a de- 
coupled system. Having determined (analytically or numerically) Z(t) from 
()14)) then we can obtain A by A 2 = dZ/dt + bZ and X^2 can be determined 
by (|15jl . For instance 

Z=fi = r 1 -l,rf = bf 1 , (ri>l), (16) 

are stationary (dZ/dt = 0) solutions (fixed points) of (fH3|) and (|TT|) for which 
reduces to the simplest (linear) harmonic oscillator d 2 F jdt 2 + i^ 2 F = 0. 
Hence A = and 

~X~ x {t) = Acos(pt), X^(t) = Asim>£) (17) 

Thus (|T7jl is a simple limit cycle solution, which is a result of a Hopf bifur- 
cation after loss of stability at r\ — 1 of the trivial solution Z =| X |= 
of (fTTf and fl~3|) . Standard stability analysis of the above fixed point shows 
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that it is stable for entire range r± > 1 so that (JT7jl is a stable cycle for all 
ri > 1. 

In the general case the Lienard equation (|T4"|) can not be solved analyti- 
cally. For its numerical integration a more appropriate form is 

d ^ = 2bZ(h-Z). (18) 

If 6 = 2/i both equations (fT4*)l and (fTSJl simplify. 

It is of interest to study the behavior of Z{t) and A 2 (t) in (|TT|) and (fT^j) 
far from the initial instant, i.e., in the posttransient period. The general 
solution of pijl with respect to Z(t) can be written as 



r . 



Z(t) = Z exp(-bt) + /*rfrexp(-6r) | X(t 

Jo 

Hence at t » t& = 1/6 (theoretically t — > oo) 

poo 

Z(t -> oo) = = / drexp(-6r) | X(t - r) 

JO 



(19) 



(20) 



The substitution of (fTHjl for Z(t) yields the integrodifferential equation 
d 



dt 



POO 

f x - / drexp(-bT) \ X(t-r) 
Jo 



X 



'21) 



valid at t —>■ oo. The integral term can be interpreted as memory of the 
process for the past. If in addition b >> 1 (in the limit b — > oo) , i.e., in the 
case of short memory, the expression (|2*U|) is reduced to 



Zoo(t) 



— ,2_ 1 ,2 



X(t) 



so that from (}2~T]) 



^ 2 - 2 ^ 2 - ¥< 



t 



oo. 



(22) 



(23) 



However (|2*3*j) is exactly the Landau equation in turbulence theory (Landau 
and Lifshitz, 1986). Its integration is an easy task. 

-i 



A\t) = A$h A 2 + (6/x - Aq) exp(-2A(t - t )) 



(24) 



A 2 , 



where A^ = [A(t )} 2 . Therefore, independently of A , A 2 (oo) = bf\ 
which means that the fixed points (fTH|l and the corresponding cycle (|17j) are 
globally stable. 
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2.2 Case of finite o (a < oo) 

Let us now assume finite a {a < oo). Then from the X, Z equations (|12jl we 
derive a new one 

§ +w =i*i'-^|i*i', (25) 

which can be rewritten as 

iWf, (26) 



or alternatively 



dt \2a 

dM 



+ 2aM = (b-2a)Z(t), (27) 



where 



dt 

M(t) = — | X | 2 -Z. (28) 
Hence in view of and ()20|) 

Moo(t) = M(t -»• oo) = ( - 1 ) r° drexp(-6r) | X(t - r) | 2 , (29) 

y2o" y Jo 

and 

= ^ I X{t) | 2 -Moo^), (30) 
la 

so that Mqo = at b = 2a. Then the first two equations (fT^jl combined with 
(|3*0)l lead to a single equation 

d^X . . dX . . 1 . l0 , . 

+ (cr + a)— + a{a - r)X + -\ X \ 2 X = 0, (31) 

valid at t —>■ oo and b = 2a. Obviously, this is a generalized autonomous 
Duffing type equation in the complex domain. In real variables it is equivalent 
to two coupled equations for Xi(t) and X 2 (t) of the same kind 

where | X | 2 = X\ + Xf , / x = n - 1, z/ = e + r 2 and X = X x + £X 2 . 

Since a — r = —j\ — iv, equation (|3*T|) has only one stationary solution 
(fixed point): X = 0. No other such solutions occur unless v = 0, known as 
"laser case" (Ning and Haken 1990). Then 

| X | 2 = 2a{n - 1), n > 1, (33) 
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and stationary oscillations in the lasers can exist. However, such constraint 
{y = 0) is not generally required for the baroclinic instability. In the latter 
case v 7^ 0. 



X(t) = Aexp(iut),uj = a ^ ,A 2 = 2a 



, e - crr 2 

(a + iy 



(34) 



is an exact limit cycle solution of the nonlinear equation (}3T]) respectively 
(JH2J), provided the term in the bracket is positive, i.e., r x > r lc = 1 + vie — 
<rr 2 )/(a + 1) 2 . 

If b 7^ 2a, the memory term aM^t) will stand in (|31j) instead of zero. 
This changes the expression for the amplitude A 2 only: 2a is replaced by b 
{b < 2a or b > 2a), thus decreasing (increasing) A 2 compared to (|3^j). 

For the case of finite a but b^Owe can use the the alternative equation 

(E3) 

^ + 2aM = -2aZ(t). 
Hence similarly to (f2T)J) and (f3T)|) 



and 



POO 

Moo(t) = -2a / dr exp(-2ar)Z oc (t - r) 



Contrary to the previous case ()2U|) the memory effect now is carried out by 
the phase variable Z(t). 

Finally, we note that equation (fH^ is identical to the last equation (JSJ) 
and equation (|3*Tj) degenerates into equation © in the real case (e = r 2 = 
0,X 2 = 0,X 1 =X,r 1 =r). 

3 The infinite Z— components version of (ED 

The paper by Booty et al. (1982) contains results for the system 

dY 00 

— = (f - 5>n)x - fl y, 

ar n=l 

^ = -6 n Z n + ^c n (X*y + XY*), (35) 



where b n , c n > 0, n = 1, 2 . . . , oo, (a, f, a) are the same parameters as in (J7J). 
However, unlike (J7|), the physical basis of the above system is the baroclinic 
instability in the geophysical fluids only. Here we are interested in (135|) from 
mathematical point of view. 
First of all if 

oo 

b n = b = const, 2J °n = c < oo, (36) 

n=l 

the system fl33|) becomes identical to ((7j) with zT = Z n . We now apply 
the approach from the section 2.2. The equation analogous to ()2(ij) here reads 

^r + fe nM n = Cn ^--lj |X|, 

where 

M n = ^ | X | 2 

Hence, at 6„ = 2<r and t — > oo we have M m — > 0, so that 

Z n (* ^ oo) = ^ | X | 2 

Further derivations follow those of section 2.2 with final result equation (p?Tj) 
with the last term on the left hand side multiplied by c = c n- Its contribu- 
tions depends on the value of c. If one adopts the model speculations from 
Booty et al. (1982) and assume c n = n~ p , p > 1, then 

oo 

C=E™~ P = C(P), P>1 (37) 

n=l 

where ((p) is the Rieman (^-function. We note that Curry (1978) and Curry 
et al. (1984) discussed system similar to (J33j) for the system ((H) 



4 Another complex Lorenz system 

Of the several complex Lorenz-like systems known in the nonlinear geody- 
namo theories (Weiss et al. 1984, Jones et al. 1985, Roberts and Glatzmaier, 
2000), the most interesting appears to be the following 

A A 

-— = \fa(a + l)DB - a A, 
dt 

dB 1 

— = iA iCA* - B, 

dt 2 

dC 

*±L = -mC - iAB, (38) 
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where unlike (J7J) all geodynamo characteristics (phase variables) A(t), B(t), 
C(t) are complex functions of the kind K(t) = Ki(t) + iK 2 (t), K = A,B, C, 
whereas D,a,m are real physical parameters and in addition a oc m oc 1. 
Thus in the real domain (f3*H|) is a sixth order nonlinear dynamical system. 
According to Weiss et al. (1984) it exhibits much richer and more fascinating 
variety of behaviour than other similar systems discussed in the cited paper. 
Below we shall concentrate our attention on it. 
Following Jones et al. (1985) we set 

A = V2X, B = (iV2/R)Y, C = 2Z/R, R = i(a + l)D/y/a. 

Then (I38j) takes the form 

^— = RX - Y - X*Z, 
dt 

— = -mZ + XY, (39) 

which closely resembles the systems (jlj and (|7jl. This fact implies that the 
subsequent results if not identical, will be similar to those corresponding to 
(JT|) and (0). For instance, at m = 2a and t » 1/m (t — > oo), the equation 
(pTTj) now reads 

^ + (a+l)^ + a(l-i?)X + i|X| 2 X = 0. (40) 

Unlike (}3T]) . X = is the only stationary solution of ()40|). On the other hand 
similar to 1)31)) and ()34j) now 

X(t) = 5 exp(^t),cu = v^f,^ = 25(D 2 - 1), (£> > 1), (41) 

is an exact periodic solution of (j4Tl|) . In details the stability and the bifurca- 
tion properties of the fixed point X = and of (|4*T]) are described in Jones 
et al. (1985). 

The physical constraints a oc m oc 1 make the previous cases (a — » oo 
and /orm^O) invalid for (}3*9"j) . 



5 Summary and conclusion 

The main results of this study concern the complex Lorenz-like systems (|7j) 
and its infinite Z-component generalization (|33|) as well as the system (|38j) . 
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arising in the laser physics and in geophysical fluid dynamics (baroclinic in- 
stability theory, geodynamo models, etc.). Some of them recover the results 
obtained by previous authors. Here we were interested mainly in the asymp- 
totic solutions of the corresponding equations: 

• At a — > oo two decoupled equations (JHJ), of known type are 
derived from the reduced system (jHJ) - (fTTj) and the solutions (JIEJ), (fTTjl 

are discussed. 

• At t — > oo a single equation (|21|) with integral (memory) term is de- 
rived. In the case of short memory (1/6 — > 0), equations (|22j) - ()24j) are 
in force. 

• At finite a (a < oo) , the generalized Duffing type equation (J3*T|) with 
exact limit cycle solution (|3*4")l is derived. Almost identical to (f3*T|) equa- 
tion (|40|) is derived for the nonlinear geodynamo model (|38|). Unlike 
(JM}, equation (jjQj) has a zero fixed point only. 

• Under the same assumptions, the generalized system (jHHjl is reduced to 
the standard form (JJJ) with Z = J2^Li Z n 

Our study open possibility for further numerical and analytical investiga- 
tion of the properties of obtained equations such as investigation of stability 
properties of attractors, bifurcation diagram and calculation of characteristic 
quantities such as generalized dimensions . This will be a subject of future 
research. 
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